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We consider a zigzag chain of coupled micropillar cavities, taking into account the polarization 
of polariton states. We show that the TE-TM splitting of photonic cavity modes yields topo¬ 
logically protected polariton edge states. During the strongly non-adiabatic process of polariton 
condensation, the Kibble-Zurek mechanism leads to a random choice of polarization, equivalent to 
dimerization of polymer chains. We show that dark-bright solitons appear as domain walls be¬ 
tween polarization domains, analogous to the Su-Schrieffer-Heeger solitons in polymers. The soliton 
density scales as a power law with respect to the quenching parameter. 


As initially shown by Kibble [T] for the expansion and 
cooling of the early Universe, and then for liquid helium 
by Zurek [5], a system undergoing a second-order phase 
transition on a finite timescale develops domains with 
independent order parameters. The Kibble-Zurek mech¬ 
anism (KZM) allows predicting the typical size of the 
domains and, therefore, the densities of the topological 
defects on their boundaries. Their scaling as a function 
of the quench rate is given by a power law, with the crit¬ 
ical exponent of the transition being determined by its 
universality class [3]. 

A very relevant system to study the KZM are the 
quantum fluids such as atomic Bose-Einstein conden¬ 
sates (BECs) formed by cooling. Indeed, quantum flu¬ 
ids support topological defects HIS], their most famous 
example being a quantum vortex, which, contrary to a 
classical vortex, like a tornado, cannot disappear ”by it¬ 
self’, via a continuous transformation. This property is 
due to the difference in vortex and ground state topolo¬ 
gies, which is guaranteed by the irrotational nature of the 
fluid described by a complex wavefunction [B] . The non- 
adiabatic cooling of such a fluid allows the development 
of topological defects obeying the KZM scaling, as con¬ 
firmed by experiments [7] and also predicted for multi- 
component BECs BM- However, solitons in ID and 
half-vortices in 2D spinor BECs maun] are only quasi- 
topological defects. Indeed, a dark soliton transforms 
into a grey and eventually disappears by simple acceler¬ 
ation, and a half-vortex can be unwound by a divergent 
magnetic field. Another system of interest to study the 
KZM are cavity exciton-polariton quantum fluids [nnm. 
Because of the finite polariton lifetime, polariton con¬ 
densation can be an out-of-equilibrium process driven by 
the condensation kinetics rather than by thermodynam¬ 
ics [m [13]. As previously pointed out[T31 [TS], the es¬ 
tablishment of a steady state by non-resonant pumping 
in an initially empty system cannot be an adiabatic pro¬ 
cess and is therefore equivalent to a quenching of the 
parameters of the system, leading to the appearance of 
topological defects. 

Another class of systems which possess topologically 


protected states are periodic lattices with topologically 
non-trivial band structures characterized by non-zero 
Chern numbers, or Zak phase, depending on their di¬ 
mensionality. The most well-known examples of such 
systems are the topological insulators [16] , Kitaev chains 
supporting topologically protected Majorana states m 
and dimer chains m- Indeed, depending on the differ¬ 
ence of the tunneling coefficients within and between the 
dimers, such chains form topologically different conduc¬ 
tion bands, characterized by a tt difference in the Zak 
phase [19]. As it was shown recently [20], the number 
of states in the conduction band depends on this phase, 
and the states which are not included in the bulk are 
localized on the edges. These edge states do not rely 
on inter-particle interactions but are topologically pro¬ 
tected: they are robust against disorder and perturba¬ 
tions. Different implementations of topologically non¬ 
trivial band structures have been studied theoretically 
and experimentally in various systems, including pho¬ 
tonics Enna, optomechanics [331 HI], excitons [3S], and 
plasmonic zigzag chains [26H29] . Optical systems offer an 
important advantage compared to the electronic ones and 
to atomic BECs, because of the facility of their fabrica¬ 
tion and the complete accessibility of the wavefunction in 
time, real and reciprocal space. Polaritonic systems were 
shaped as molecules and lattices [3(1132] . Schemes for 
creating polariton topological insulators have been pro¬ 
posed [53H55] . While KZM, as a universal mechanism, 
has already been widely studied in various systems shar¬ 
ing some common properties with our proposal, such as 
zigzag ionic chains [371140] . where the phase transition 
and the topology correspond to physical arrangement of 
atoms, none of these possess the same key ingredients. A 
topologically non-trivial polaritonic chain therefore ap¬ 
pears as an ideal system to study the complex interplay 
of topological ordering and KZM m- 

In this work, we describe polariton BEC in a zigzag 
chain of polariton micropillars with photonic spin-orbit 
coupling (SOC) [331143| 133]. As a result, the polariton 
band is characterized by a non-zero Zak phase and the 
chain supports topologically protected edge states. These 
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results are also valid for Rashba SOC present in atomic 
condensates [44]. We show that with a focused non¬ 
resonant excitation spot condensation occurs on the edge 
states with polarization determined by the Zak phase. 
When the system is excited homogeneously, a gas of dark- 
bright solitons (a spinor version of Su-Schrieffer-Heeger 
(SSH) solitons) is formed via the KZM. Here, the TE- 
TM SOC is crucial, allowing a homogeneous condensate, 
contrary to fc-linear Rashba SOC. We demonstrate that 
the soliton density follows a power law with respect to 
the quenching parameter - pumping intensity. The scal¬ 
ing exponent for different system parameters is close to 
1/4, in agreement with the mean-field KZM theory for a 
ID system. 

Existence of topological edge states. We first consider 
a zigzag chain of coupled OD modes neglecting the spin, 
as usually done in theoretical analysis of the electronic 
dimer chains. Let us call the first pillar in the chain ”a”, 
and let it be ’’below” the second pillar ”b”, so that the 
first link is oriented at 45° (see Fig. 1), which we will call 
’’diagonal” direction (D), while the perpendicular direc¬ 
tion shall be ’’anti-diagonal” (A, 135°). The pair ”ab” 
forms the unit cell. Following the definitions established 
in the previous works |20j , the tunneling constant in the 
first link (within the cell) is called t', while the tunneling 
in the second link (between cells) is t. The corresponding 
tight-binding Hamiltonian reads (m is the cell number): 

H = '^t'b'l^dm + tdl^_^_lbm +H.c. (1) 


where a and b operators act on the corresponding pillars. 

Let us now consider that these OD modes are con¬ 
stituted by photonic micro-pillars obtained by etching 
a planar cavity [IS]. Each pillar ground state has two 
polarizations, which we assume to be degenerate. On 
the other hand, the optical eigenmodes of the cavity are 
TE and TM polarized and have different effective masses 
pS] . This makes the tunneling coefficients polarization- 
dependent [Ml SH133] and different for the polarizations 
oriented longitudinal and transverse with respect to the 
link. We therefore have t < t' for D-polarization, for 
which the first link (labeled by t') is longitudinal (Fig. 
l(a,d)), and t > t' for the A-polarization, for which the 
same link is transverse (Fig. l(b,e)). The relative differ¬ 
ence in the longitudinal and transverse tunneling coef¬ 
ficients for typical parameters of a polariton micropillar 
lattice can be of the order of 10% [34]. This difference of 
the tunneling coefficients is equivalent to the dimerization 
of polymer chains, but associated with the polarization 
of the states. The corresponding dimers are shown with 
black dashed lines in Fig. 1. 

For an even number of pillars, one can directly apply 
the SSH theory developed for polymer chains [15]. The 
dispersion of such system contains two conduction bands, 
below and above the single-pillar energy, chosen as the 
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FIG. 1: (color online) A scheme of a zigzag chain with even 
and odd number of polariton pillars: ”a” and ”b” form a unit 
cell marked in dark grey. a,d) Diagonal polarization: t' > t, 
b,e) Anti-diagonal polarization: t' < t. c,f) Energy band of 
30(c) and 31(f) pillar chains obtained from the tight bind¬ 
ing Hamiltonian. The color shows the diagonal polarization 
degree of the states. 


zero reference. The existence of the edge state in this 
case is determined by the Zak phase, which is an analog 
of the Berry phase [37] defined on a unit cell of a size d 
for a Wannier function Unk (a:), integrated over a given 
band n [T9] . 


Cn = 




{x)i ^Unkix) d^dk 

ok 


( 2 ) 


The Zak phase is determined by the ratio of the tun¬ 
neling coefficients within and between the dimers [20] : 
Cn = 0 if t'/t > 1, and = tJ", if t'/t < 1. The topo¬ 
logical transition = 0 O Cn = occurs a.t t = t'. If 
the Zak phase is tt, the number of states in the bulk is 
less than the number of pillars: N = M — 2, and the re¬ 
maining states, whose energy is that of uncoupled modes, 
are localized on the edges of the chain. If the Zak phase 
(n = 0, the number of states in the bulk is equal to the 
number of pillars JV = M, and no edge states appear. 
The advantage of the optical systems is that the Zak 
phase can be measured directly m- 

In our system, as can be deduced from Fig. 1(a) and 
(b), a pair of edge states does exist in one polarization 
(anti-diagonal for our parameters), and does not exist 
in the other. The result of the diagonalization of the 
Hamiltonian (1) with spin (see [35]) for a chain of 30 pil¬ 
lars {5J — 0.3J for visibility) is shown in Fig. 1(c), with 
polarization shown with color. We see that the polar¬ 
ization states are interleaved, the lower band ends with 
polarization D, and the upper band begins with D, so 
the edge states (seen in the gap) are both necessarily A- 
polarized. Rashba SOC gives a similar result. For odd 
number of pillars, if t/F ^ 1, one can redefine a dimer so 
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that the Zak phase will be tt and a state will appear on 
the edge which contains a unit cell broken by the bound¬ 
ary. For polaritons, an important consequence is that for 
one polarization (D), the edge state is on the right edge 
of the chain (Fig. 1(d)), and for the other polarization 
(A) the edge state is on the left edge of the chain (Fig. 
1(e)). Calculation yields Fig. 1(f), where for all states, 
including the edge ones, the polarization is interleaved. 

Overall, whatever the number of pillars in a finite 
zigzag chain, because of the polariton SOC there are al¬ 
ways two edge states in the system, either having the 
same polarization when the number of pillars is even, or 
being cross-polarized when the number of pillars is odd. 
The edge state polarization is always orthogonal to the 
axis linking the two last pillars at the chain edge. 



FIG. 2: (color online) Calculated spatial images of the differ¬ 
ence in diagonal polarization emission lo— I a'- a) edge states 
in a chain with an even number of pillars; b) edge states in a 
chain with an odd number of pillars; c) Emission of the con¬ 
densed states under localized pumping (marked P). Opposite 
diagonal polarization is observed on opposite sides. 


We confirm the predictions of the analytical tight- 
binding model by solving numerically the spinor 
Schrodinger equation on a grid to find the eigenstates: 

B,Ji = -AAfc + ^(|-=F.|-) + (3) 

where = (i/;+(r,t),^/i_(r,t))^ are the two 

circular components of the wave function, j3 = 
nil ^ while mi t are the effective masses 

of TM and TE polarized particles respectively and m = 
2 {nit — ’mi) /rntmi', mt = 5 x 10“^mo, mi = 0.95mt] toq 
is the free electron mass; U{r) is the potential of the pil¬ 
lars describing the confinement of polaritons in the chain 
beyond the tight-binding model [M] . The results of these 
calculations are presented in Fig. 2, showing the spatial 
images of the difference between the diagonal polariza¬ 


tions Id — Ia for both even (a) and odd (b) number of 
pillars. The localization length is discussed in supple¬ 
mental material |48j . 

Condensation on localized edge states A very effective 
way to excite these localized edge states is to create a 
polariton condensate using focused non-resonant pump¬ 
ing. This technique allows creation of strongly out-of- 
equilibrium states, typically the states showing the best 
spatial overlap with the localised excitonic reservoir in¬ 
duced by the pump [32l |49l HD] . In lattices, the repulsive 
potential induced by the excitonic reservoir becomes at¬ 
tractive for particles with negative effective mass at the 
band edges, which leads to the condensation on localized 
gap states bound to the reservoir EaEOj. One expected 
peculiarity of the zigzag chain with a local pump is that 
the polarization of the localized mode where the conden¬ 
sation occurs is entirely fixed by the chain topology and 
the position of the pump, and does not rely on a symme¬ 
try breaking process. To demonstrate this predicted fea¬ 
ture, we model polariton condensation using the Hybrid 
Boltzmann-Gross Pitaevskii equation which includes re¬ 
laxation mechanisms EllillliSI- For a thermal excitonic 
reservoir, the model can be reduced to: 

= V-T (4) 

+U'ip± - llV’i + ((ffii + *7(»^))^± +X)exp ) 

where the parameters (values from ref. |32j 1 other than in 
Eq.(3) are: A - the kinetic energy relaxation term, Ur - 
the reservoir potential amplitude, a - the reservoir width, 

T the lifetime, x ~ fh® Gaussian noise term included to 
describe the spontaneous scattering 7 ( 71 ) is the 

saturated stimulated scattering rate from the reservoir, n 
is the total polariton density. We neglect the interactions 
within the condensate here. A pump located close to the 
edge will excite the unique localized mode at the edge of 
the chain. If the pumping spot is located in the bulk, the 
potential of the reservoir cuts the lattice into two smaller 
chains, and the same reasoning as above applies to each 
of them, leading to the condensation at their respective 
edge states. The results of the simulations for the pump¬ 
ing spot located in the middle of a chain (which allows 
to check all predictions simultaneously) are presented in 
Fig. 2(c), showing the difference between the intensities 
of the diagonal polarizations Id — IA of the light emitted 
from the system above the condensation threshold. The 
condensation indeed occurs on the localized edge states 
on both sides of the spot, with polarization controlled by 
the condition on the Zak phase = tt. 

Spontaneous formation of dark-bright solitons via the 
Kibble-Zurek mechanism The consequences of the non¬ 
trivial topology of the system are truly revealed under 
homogeneous non-resonant pumping. While the previ¬ 
ous results could be verified by purely linear measure¬ 
ments in this part we study the condensation of 
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FIG. 3: (color online) a,b) Difference in diagonal polarization 
intensities for 2 examples of polarization textures in a zigzag 
chain. Colors as in Fig. 2. c) Average number of solitons as 
a function of effective pumping fitted by power law. 


polaritons, accompanied with the emergence of a spinor 
order parameter. As any second-order phase transition 
occurring on a finite timescale, it is described in terms of 
parameter quenching, responsible for the KZM-type for¬ 
mation of topological defects HHH]. In a spinor system, 
the large phase fluctuations at the early stage of a con¬ 
densation process are associated with spatial fluctuations 
of the polarization, governed by the phase difference be¬ 
tween the spin components. The polarization domains 
correspond to the dimerization domains in the SSH pic¬ 
ture, and the domain wall between them is equivalent 
to the SSH soliton [T5]. These dark-bright solitons sep¬ 
arate two polarization domains, both equally stable and 
characterized by a fixed difference tt in their Zak phase. 
The stability of a dark-bright soliton, which cannot be 
destroyed by acceleration, is its most important feature: 
once formed, it does not disappear, contrary to a scalar 
dark or grey soliton. This stability requires the interac¬ 
tions in the condensate to be smaller than the TE-TM 
splitting [48]. It allows to detect such objects in cw exper¬ 
iments, being advantageous with respect to the previous 
proposals of the KZM studies with polaritons miiii, 
requiring single-pulse experiments. 

Fig.3 shows the results of simulations based on the nu¬ 
merical solution of the equation (4) with a homogeneous 
reservoir potential (ct = oo) and interactions |3H]. Panels 
(a,b) show the difference between the intensities of the 
diagonal polarizations Id — Ia oi the light emitted by a 
condensate formed under weak and strong pumping re¬ 
spectively. In panel (a), two domains A- and D-polarized 
corresponding to a single dark-bright soliton are visible. 
Panel (b) shows 6 polarization domains and 5 domain 
walls (for a movie see supplementary |48j ). In our numer¬ 


ical experiment, we do not change the temperature of the 
system, as in the classical KZM, but rather turn on the 
pumping and fill the system with particles, changing the 
critical condensation temperature, but keeping the sys¬ 
tem temperature (controlled by y and A) constant. The 
quenching time is controlled by the pumping intensity 
Tq^ (X P. This scheme, while being simpler and ubiqui¬ 
tously present in all polariton experiments, fits the KZM 
scheme, because the relative temperature e = (T—T^j/Tc 
at threshold changes linearly with time (see [48| ). 

Indeed, the Gaussian noise y is uncorrelated in time, 
and so creates a frequency-independent population of 
particles. However, the energy relaxation term propor¬ 
tional to A describes energy-dependent decay acting on 
these particles. The resulting spectral density \tp{E)\'^ 
for the polariton state of energy E can be obtained as: 
\ip{E)\'^ oc y/r, where T is the total decay rate, com¬ 
posed of energy-independent Tq (ground state lifetime) 
and energy-dependent relaxation Ta = AE [51], giving: 


\'ip{E)\'^ oc 


X 

To+AE 



(5) 


The linear part of the spectrum at low energies can be 
interpreted as a Boltzmann distribution function with 
an effective temperature T = Tq/A. Varying the relax¬ 
ation efficiency, we can change this effective tempera¬ 
ture: the better is the relaxation, the lower is the tem¬ 
perature. According to KZM, the average density of the 
topological defects Ugoi in a BEG scales as the inverse 
healing length kTi where z/ is a scaling ex¬ 

ponent. In the mean-field approximation, zz = 1/2 (be¬ 
cause ^ = H/'/2anm, where a is the interaction con¬ 
stant), and the dynamical exponent is z = 2 [3] [S] 
for relaxation linear in energy, which allows to write 
nsoi = oc Thus, Usoi oc 

oc 7^/'*: a scaling exponent of 1/4 is expected for 
the mean-field universality class. Fig. 3(c) shows the 
number of solitons appearing in a chain of 40 pillars ver¬ 
sus the effective pumping intensity (7 — re//)/ro, where 
Te// « S.STq is the effective decay rate accounting for 
A. Each point is obtained as an average of 10 simula¬ 
tions. The power law fit is compatible with the scaling 
exponent 1/4 expected for KZM. 

To conclude, we have demonstrated that the Zak phase 
plays a crucial role for the description of condensation 
in ID zigzag chains of polariton pillars. Because of the 
SOC, such chains always exhibit exponentially localized 
edge states. Similar results can be obtained for Rashba 
or Dresselhaus SOC in atomic BECs. Under homoge¬ 
neous pumping, dark-bright solitons appear between the 
domains of orthogonal polarization via the Kibble-Zurek 
mechanism. We extract numerically the dependence of 
the soliton density against the quenching parameters 
and find it in agreement with the analytical predictions. 
These domain walls can also be created by quasi-resonant 

















5 


excitation, for example, using Gauss-Laguerre beams fo¬ 
cused on a pillar chain (above the bistability threshold) 
m- They can also be manipulated using the electrically 
controlled in-plane effective magnetic fields [SS], which 
might allow to design optical race-track memories m- 
We acknowledge discussions with M. Glazov, A. Amo, 
and J. Bloch and the support of ITN INDEX (289968) 
and ANR Labex GANEX (207681). 


SUPPLEMENTAL MATERIAL 

Tight-binding Hamiltonian with polarization 

The matrix of the tunneling couplings in a single zigzag 
chain in the tight-binding approximation taking into ac¬ 
count the spin-dependent tunneling is given by 


Fk 


{ hJ f^SJ\ 
Kf^sj f^j ) 


( 6 ) 


Eig. S3, where the numerical simulations (using the 2D 
model described in the main text) are compared with the 
analytical expression obtained in the framework of the 
tight-binding approximation. The good quality of the 
approximation in this case is due to the direct tunneling 
of photons within the confinement potential. 



where complex coefficients /k,/^ are defined by: 

/k = exp(-ikd,^), = exp(-i [kd,^ =F 2^?]), 

and if = ±7r/4 is the angle between horizontal axis and 
the direction to the jth nearest neigbor. Each pillar in 
a chain has only 2 neighbours, contrary to the planar 
graphene case. J is the polarisation independent tunnel¬ 
ing coefficient, whereas 5J \s the SOI-induced polariza¬ 
tion dependent term. The results of the diagonalization 
of the Hamiltonian are shown in Eig. l(c,f) of the main 
text, where the edge states appear in between the two 
dispersion bands. The polarization degree of the bulk 
states on this figure is smaller in the case of odd num¬ 
ber of pillars, because the extra pillar affects the overall 
equivalence of the two diagonal polarizations. 

Localization length of the edge states 

The inverse localization length k can be estimated by 
requiring that the corresponding wavefunction vanishes 
on the opposite edge: 

t' sinh/t (M + l)a = t sinh uMa 

where a is the length of the unit cell, and M is the 
number of pillars in the chain. In the limit of suffi¬ 
ciently large difference of t and t', the equation reduces 
to t'/t « exp (—Ko) which means physically that we deal 
with a particle whose mass is given by the larger of the 
coefficients t, which decays inside a barrier of a height t' 
(which is the energy gap between the two bands). 

The numerical simulations also confirm that the decay 
of the edge state is indeed exponential, as can be seen in 


FIG. 4: (color online) The exponential decay of the edge state; 
solid line - numerical simulation, dashed line - analytical so¬ 
lution. 

The exponential localization of the edge states plays 
a very strong role for the condensation on these states 
in case of localized pumping, as explained in the main 
text. In case of pumping on the edge, the localized edge 
states, which always exist in polariton zigzag chain, as 
shown above, will have a much larger overlap with the 
reservoir than any of the propagative bulk states, with 
the overlap ratio being L/l^ where L is the length of the 
chain and I = k~^ is the localization length of the edge 
state, discussed above. 

Formation of a condensate and the Kibble-Zurek 
mechanism 

In this section we discuss the simulations and the corre¬ 
sponding experiments required for the study of the topo¬ 
logical defect formation via the Kibble-Zurek mechanism 
under non-resonant pumping. The system is composed 
of a homogeneously pumped exciton reservoir and po¬ 
lariton states. As written in the main text, we describe 
the system by the spinor Gross-Pitaevskii equation with 
phenomenological terms describing decay, energy relax¬ 
ation, stimulated scattering from the reservoir into the 
condensate, and the noise, describing the spontaneous 
scattering: 

*^^ = -(1-*A)£aV^±+/3(J|t*z|)V (7) 

+a\tp±f iIj± -f UiIj± - ^tp± + {Ur + ij{n)) + x 




6 


Here, a is the interaction constant for polaritons hav¬ 
ing the same spin. For opposite spins, the interactions 
pass through the dark exciton states, and the correspond¬ 
ing constant is negligibly small, as a second-order correc¬ 
tion. Thus, the interactions are strongly spin-anisotropic, 
which allows to consider the two spin components as be¬ 
ing almost decoupled during their formation. The stim¬ 
ulated scattering term +i'y{n)ij} in the Gross-Pitaevskii 
equation leads to the exponential increase of the par¬ 
ticle density with time, until the saturation is reached. 
The saturation is obtained by the inclusion of exponen¬ 
tial saturation factor 7 (n) = 70 exp(—It is 
the constant prefactor 70 which determines the speed 
with which the condensation threshold is passed e{t) = 
(T — Tc)/Tc{t), and therefore, the number of topological 
defects via the Kibble-Zurek mechanism. In this type of 
experiment, it is not the system’s effective temperature 
T, but the critical temperature Tc which changes with 
time, while the system passes the condensation thresh¬ 
old. In ID, Tc oc n^, where n oc exp( 7 t), which, after the 
linearization of the exponent, gives at threshold a linear 
behavior e{t) = t/rg, where tq oc A/ 2 ro 7 , allowing to 
test the scaling behavior as a function of 7 . 

In experiments, it is the pumping intensity which con¬ 
trols the reservoir density, which, in turn, controls the 
stimulated scattering efficiency. The dominant relaxation 
mechanism is usually the exciton-exciton scattering, and 
thus the probability for the reservoir particle to scatter 
into the condensate is proportional to Wo{l+Ncond)N!^es^ 
where Wq is a constant prefactor proportional to the 
exciton-exciton interaction strength. Varying the reser¬ 
voir density by changing the pumping, one can therefore 
change the stimulated scattering rate 70 oc in exper¬ 
iments. Of course, this also changes the saturation value 
for the condensate density, but the Kibble-Zurek mech¬ 
anism depends only on the parameters important at the 
transition point. We have specifically checked in numer¬ 
ical simulations, where the parameters can be changed 
independently, that the saturation density Umax does not 
affect this mechanism, if the density remains sufficiently 
small for the TE-TM coupling to remain dominant over 
the interactions. If the interactions dominate the spin- 
orbit coupling, the topological gap closes, the dark-bright 
solitons are replaced by less stable half-solitons, and, fi¬ 
nally, a homogeneous linearly polarized condensate, cor¬ 
responding to the ground state of the system, emerges. 
In our simulations such behavior occurs if the pumping 
intensity is increased significantly above the range shown 
in Fig. 3 and Fig. S2. 

The behavior of A as a function of the reservoir density 
also has to be taken into account to describe the exper¬ 
iments properly. If a linear dependence is assumed for 
exciton-exciton relaxation mechanism, Tq^ oc 7 /A oc T, 
where P is the pumping of the reservoir, which means 
that dependence on P in experiments can be expected to 
be the same as the dependence on 7 in numerical simu¬ 


lations. 


Additional results on KZM 

In this section we present additional results with KZM 
simulations. We have compared the formation of topo¬ 
logical defects for two different values of the relaxation 
constant A, determining the quench time tq oc A. Fig. 
S2 shows that, as expected, a faster quench (smaller A) 
leads to a higher number of solitons, but the power law 
remains the same, with a scaling exponent 1/4. A smaller 
system size of 19 pillars, as compared to Fig. 3c) of the 
main text (40 pillars), was used for these simulations, in 
order to shorten the calculation time. 



(Y-re«)/r„ 


FIG. 5: (color online) Average number of solitons as a func¬ 
tion of effective pumping for different relaxation efficiency: 
A = 0.06 (black) and A = 0.03 (red). Power law fits are 
shown with solid lines. 

We note that both in the main text and in the sup¬ 
plemental material, the fit gives a slightly higher scaling 
exponent of about 0.28 instead of 0.25. Current precision 
of the simulations does not allow to distinguish between 
the two values, but future studies might show if this de¬ 
viation from the analytical predictions is significant, and 
if it is due to the particular (for example, topological) 
properties of the system. 

Finite-size effects take place in our system for higher 
values of pumping, when the distance between the soli¬ 
tons becomes comparable with the period of the chain. 
In this case, topological protection is reduced, and soli¬ 
tons recombine, which leads to their average number be¬ 
ing much smaller than the one predicted by scaling law. 
These effects begin approximately at y/Fp = 20 for 19- 
pillar chain and at 7 /F 0 = 35 for a 40-pillar chain. This 
is why it is difficult to obtain a large number of solitons 
in the chain. 
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VIDEO 

The supplemental Video file demonstrates the polari- 
ton condensation in a zigzag chain of micropillars. The 
colormap corresponds to the difference between the in¬ 
tensities of the diagonal polarizations, (Id — J^), as in 
the figures 2 and 3 of the main text. One can observe 
the formation of several polarization domains and the 
establishment of a steady state. 
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